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0. Introduction 

One of the main themes in discrete geometry has always been the study of op¬ 
timally dense packings of bodies, in regions of finite and, especially, infinite volume 
- the latter being the situation with which we will be mainly concerned. There 
is a large literature for packing in Euclidean spaces E n . but much less for packing 
in hyperbolic spaces EI n ; see [GrW] and references therein, as well as the classics 
[Fe5, Rog]. The difficulties in hyperbolic spaces have been documented in many 
papers dating back at least to the early 50’s (see [Fel-5], [Bol-2], [BoF], [FeK], 
[FKK], [Kup]), and are well understood to be related to the phenomenon of the 
exponential rate of growth of the volume of a sphere with respect to radius. 

In general it is difficult, even in E" . to actually determine the optimal packing 
density even for simple shapes; for instance this is unknown even for spheres for 
dimension n > 4. For certain polyhedra such optima are occasionally computable, 
and one outgrowth of this, since the late 60’s, has been the subject sometimes re¬ 
ferred to as “aperiodic tiling”; see [Rad] and references therein. That subject is 
concerned with the geometric symmetries of those packings (tilings) which achieve 
optimal density, symmetries of unusual form, which are studied through probabilis¬ 
tic techniques. 

In hyperbolic space the difficulties in analyzing optimal density are, as hinted 
above, much worse; not only is it hard to determine the optimal density, say for 
spheres, it has even been hard to decide whether this is a meaningful question. 

Finding a meaningful, and computationally useful, way to analyze optimal 
density in HP was the subject of [BoR, Bow], in which some of the probabilistic 
methods of aperiodic tiling were employed. The ergodic theory approach to den¬ 
sity in [BoR], summarized below, is a modification of an approach based on Mass 
Transport suggested to us by Oded Schramm (see [BeS]). We found, using ergodic 
theorems of Nevo et al, that the quantity he suggested as a “density” is in fact, in 
a statistical sense discussed below, the highest true density - the limiting fraction 
of space covered by bodies in a fixed packing. The papers [BoR, Bow] did not, 
however, convincingly address the question of optimally dense packings themselves; 
that is, a formalism for analyzing the density was found, but not a way to address 
the packings that achieve that optimal density. That is one of the main aims of this 
paper. 

The other major objective of this paper is to present a series of examples, some 
new and some which are variants of historically important examples, adjusted to 
reveal significant connections. Although one source of difficulty in the subject is 
well understood, the rate of growth of volume with diameter, these examples point 
to a very different source, related to structures in the space of packings. 

I. Some troublesome packings 

We will be analyzing the density of certain subsets of Euclidean n-dimensional 
space E n or hyperbolic n-dimensional space HP of curvature — 1; we let § stand for 
any of these spaces. The subsets of § whose density we will consider will generally 
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be “packings” of (infinitely many) “bodies” /3j, where a packing is a collection of 
bodies with pairwise disjoint interiors, and a body is a compact, connected set which 
is the closure of its interior. One of the features of our analysis will be an emphasis 
on distinguishing between the density of the packing versus the density of the set 
which is the union of the bodies in the packing; that is, it will be significant to 
maintain the individuality of each of those bodies. 

Our main focus will be on the “densest” packings possible by the given bodies, 
and this requires examination of the primitive notion of density. If we were packing 
a region S of finite volume by the bodies (3j , the density of such a packing would be 
unambiguous - the fraction of the volume of S covered by the bodies - but density 
must be defined more subtly for packings of a region, such as S, of infinite volume. 
The most widely accepted [FeK] primitive notion is that the density of a packing 
V of § should be obtainable by choosing a family of finite volume regions Sk, with 
Sk C Sk+i and U k.Sk — §, and the density of V should be 


lim 


vol{V fl Sk) 


( 1 ) 


vol(Sk) 

where vol(-) denotes volume in § and V fl Sk denotes the portion of Sk covered by 
bodies in V. We would want the density to be reasonably independent of the family 

S k - 


It is worth noting that the limit in (1) can easily fail to exist. Consider the 
sequence {Dj j > 1} of closed disks in E 2 , Dj of radius 2 :l and centered at the 
origin. Let Pj be the annulus Dj/Dj- 1 between successive disks, and let S be 
the union of those Pj with j > 2 even. If we try to define the density of S using 
the expanding regions Sk = Dk, the sequence of local or approximate densities 
vol(S fl Sk)/vol(Sk) would not have a limit as k —> oo, due to oscillation. (We 
could get the same qualitative result by replacing our region S by its intersection 
with some simple packing of disks, such as the packing of unit diameter disks whose 
centers have integer coordinates.) 

Even though there are packings without a well defined density there is no real 
difficulty in defining optimal density of packings in Euclidean space. In fact we now 
show how to construct densest packings of Euclidean space. Let § = E" and let Sk 
be a cube centered at the origin, with edges of length k aligned with the axes. For 
any k > 0, let Vk be a packing by (congruent copies of the bodies in) B = {Pj} 
such that all bodies in Vk intersect Sk and vol(Vk fl Sk) is optimally large. (Such a 
packing is easily shown to exist by a simple compactness argument [GrS; p. 154].) 
For any packing V in § we define 


MV) = vol{ %\ p) 
vol(Sk) 

(2) 

dk — ma xdk(V) 

'p 

(3) 


d = lim sup dk ■ 

k —>oo 


(4) 
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(Hindoo d k exists but we do not need this fact.) 

At this point it is convenient to have a space Eg of all possible packings of § 
by the bodies (3 3 , equipped with a metric topology such that a sequence of packings 
converges if and only if it converges uniformly on compact subsets of §. We will 
spell this out in section II b, but assume for now such a space makes sense and is 
in fact compact. Then we let Voo be an accumulation point of {V k }• 

The following is a simple observation. 

Lemma 1. dk(gVoo) —> d as k —> oo for every fixed rigid motion g. 

Proof. The main estimates needed are the simple facts, for k' > k: 

voi{v k n s k ) > voi{Vk> n s k ) (5) 

voi(v k , n s k ) > voi(v k n s k ) - [ k n -(k- c) n ], (6) 

where C is larger than the diameter of any body in B. Equation (6) holds because 
if it did not one could arrive at a contradiction by altering V k ' as follows. First 
replace the bodies of V k > that are completely contained in S k by the bodies of V k 
that do not overlap the other bodies of P k > (i.e. that do not overlap any body of P k > 
that overlaps the complement of S k ). Note that the volume of bodies of P k that we 
have introduced is at least as large as the right hand side of (6). Since vol(V k >) is 
as large as possible, this operation could not have increased its volume. This proves 
(6). Since Poo is a limit of V k we get that (6) holds if P k > is replaced by Voo- 

Finally, if k rn is a sequence such that d km —> d as m —> oo: 

\d - d km {gV oo)| < | d km - d\ + \d km - d km {gV oo)| . 

= I d km -d\ + | d km (V km ) - dkJgVoo)] 
and | d km (Vkm) ~ d km {gV, oo)| -> 0 as m -> oo from (6). | 

Thus, in Euclidean space optimally dense packings V exist for any collection B 
in the sense that their density defined by (1) exists, and is as large as that for any 
packing. 

As we shall see, the above technique does not extend to § = HP and therefore 
some other method must be used to define optimal density in HP. Before exhibiting 
such a method, we present some examples to highlight some differences between 
hyperbolic and Euclidean packings. 

Example 1 (half-space). 

Consider the half space region S, defined, in the upper half plane model of the 
hyperbolic plane, as the set of points (x,y) with x > 0. If we try to define the 
density of this region by circles all expanding about a common center c, it is easy 
to see that the density would depend on c, with any value strictly between 0 and 1 
being obtainable for appropriate c. This suggests that we will want the origin, used 
for the expanding regions in (1), to be arbitrary. 
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Example 2 (stripe model). 

We now give a simple example of a region S in the hyperbolic plane such that, 
when we try to define the density of S relative to a sequence of circles expanding 
about some point, we get the kind of oscillation we found in the Euclidean annulus 
example. We define the “stripe model” in the (upper half plane model of the) 
hyperbolic plane, where the stripes are the regions separated by the horocycles hj , 
j e Z, defined by y — y 3 = e^ +1 / 2 ^ w , where fixed W » 1 is to be specified. 
These curves are equidistant by W in the hyperbolic metric. We call those stripes 
separated by li 2 3 and /i 2 j+i “black”, and the others “white”, and we declare the 
region S of interest to be the union of the black stripes. 

Consider the circle with hyperbolic center c = (0,1) and hyperbolic radius 
R — (N + 1/2) W, where N » 1 is to be specified. We will use the following 
relations between the hyperbolic center (if, K ) and hyperbolic radius R of a, given 
circle and its Euclidean center (h, k ) and Euclidean radius r: 


h = H, k 2 — r 2 = fi 2 , r = ktanh(R). (8) 


So our circle has Euclidean center (0, cosh[f?,]) and Euclidean radius sinh(f?). 

We will show that, if N is even, the area inside the circle, of the black stripes 
is larger than that of the white stripes; in particular, each black stripe, between hj 
and hj+i, j < N — 3, is larger (by a factor 2) than that of the neighboring white 
stripe above it (between hj+ 1 and hj+ 2 ), and therefore the area of the circle is at 
least 2/3 black. 

For —N — 1 < j < N — 1, the area Aj of the stripe between hj and h J + \ is: 




rVj+i /* [2j/ cosh(ii) —1—j/ 2 ] ^ ^ 


'Vi 


' i — dx dy 

— [2y cosh(H) — 1 — y 2 ) ^ U 


ryj+i 2[2y cosh(R) — 1 — y 2 ] 2 


(9) 


'Vi 


y“ 


dy. 


For —N < j < N — 2 the leading behavior as N, W —> 00 (and recalling that 
R = [N + 1/2]W), is 


Aj ~ 


" y 3 + 1 2 y 1 / 2 e R / 2 


dy 


'Vi 


( 10 ) 


4e [R/2-(j + l/2)W/2] 


where a ~ b means-» 1 as N, W —> 00 . So 

b 


A 


A 


J ~ e w / 2 . 


3 +1 


(ii) 


For j = —N — 1 we have: 
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so 


A-n-i 


r— 2 ( e R y-l^ 

Je-R V 



(z- l ) 1 / 2 


dz 



—7 dz 

z 2 


> 


e 


R 


( 12 ) 


A-n-i 

A-n 


> i P w /2 

~4 ' 


(13) 


Finally we note that | e w / 2 can be made as large as desired, in particular larger 
than 2, which completes the argument that the relative densities of the set S of 
black stripes does not have a well defined limit. 


The example of the stripe model in the hyperbolic plane, where the stripes are 
all of equal “width”, is more unsettling than the example of annuli in Euclidean 
space discussed above, where in a sense the oscillation was more obviously built in. 
We will see below that this stripe model is only a simple version of a well known 
disk packing. 

There has been another common way to compute or estimate the density of 
packings in Euclidean spaces, using tilings associated with the packings, and the 
relative densities of the bodies in the tiles. (A tile is a homeomorphic image of the 
closed unit ball, and a tiling is a packing by tiles for which the union of the tiles is 
the full space §.) We emphasize that this is an attempt to reduce the intuitive global 
idea of density, which involves taking a limit of approximate densities in expanding 
regions of finite volume, to a more local notion. As a significant example of this 
approach we note an elegant proof [Fel, Rog] of the optimal density for packings of 
equal disks in the Euclidean plane. The proof uses the Vorononi cells of the bodies 
of a packing, where the cell for a body j3 is the set of all points p 6 § as close 
to (3 as to any body of the packing. The proof shows that the relative density in 
its Voronoi cell of any disk of any packing is bounded above by that of any of the 
Voronoi cells in the obvious hexagonal packing. This argument was extended to 
sphere packings in 8 by K. Boroczky, who showed [Bo2] that the relative density of 
any sphere of any packing of § in its Voronoi cell is bounded above by the relative 
density associated with that of a regular simplex. (See [FeK] for details.) Such 
relative densities in tiles of associated tilings have remained an important tool in 
analyzing optimal densities of sphere packings in Euclidean spaces [FeK, Bez], 

Example 3 (tight radius packings) 

In hyperbolic space, particularly the plane H 2 , the above method of estimating 
or computing a density of sphere packings through an associated tiling has been 
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used convincingly for the special case of disks of “tight” radius. The radius r of a 
sphere in S is called tight if the regular simplex of side length 2r admits a (full-face 
to full-face) tiling of 8. In H 2 this is the case if and only if the equilateral triangle 
of edge length 2r has angles of the form 2ix/n for some m > 7, in which case 

2r = 2 r m — cosh -1 [cot ( — ) cot( —)] (14) 

m m 

and clearly r m —> oo as m —> oo. For disks with tight radius r m the obvious 
“periodic” packing, in which each disk is surrounded by m disks touching it, has 
a well defined density in the sense that, besides the method using Voronoi tilings, 
any reasonable way to compute the density would give the same value (namely 
[3csc(7r/m) — 6]/[m — 6] [Fe5]), in particular any limit of the form (1) [BoR], 

Example 4 (Boroczky’s packing). 

There is an influential example due to Boroczky [Bol] which points out a diffi¬ 
culty in using relative density in tiles to define the density of at least some packings 
in hyperbolic space, even some which are rather symmetric. Place disks in the up¬ 
per half plane model of the hyperbolic plane with Euclidean centers at those points 
with coordinates 


{(e 2j+ Hk + 1/2), e 2j+ ^)\ j,keZ}. (15) 

The connection between this and the stripe model is simple: we are placing the 
disks equally spaced in the black stripes (and we are taking the value W = 1 for 
the width of the stripes). See Figure 1 for a picture of the packing, which includes 
some horocycles and geodesics to help understand the structure. 

In Figure 2 we see the same packing with two congruent tiles in dark outline. 
For each tile consider the tiling of the plane made by congruent copies of the tile, as 
follows. First produce copies of the tile by the congruences: (x,y) —> (x + mw,y ), 
m E Z, where w is the Euclidean width of the body. This fills out a black and 
white stripe. Then produce, from these, more copies of the tile by the congruences: 
(x,y) — > (e 2m x, e 2m y), m E Z. Together these copies of the original tile will cover 
the whole plane. The two tilings made this way, one from each of the tiles in Figure 
2, are both simply related to the same packing of disks. The punchline is, the 
tiling made by starting with the tile on the left in Figure 2 would suggest assigning 
a “density” of the packing of disks twice the value suggested by the tiling made 
by starting with the tile on the right! We repeat the point that using a tiling to 
compute the density of some packing, thus making the computation more local, is 
useful in Euclidean spaces but is less convincing in hyperbolic spaces. 

We now return to the question of a definition of optimally dense packings of 
HP. As we say above, for packings of Euclidean space the notion of densest packings 
is easy to clarify, and one way to understand this is through the computation of the 
ratio /(p, a) of volumes of concentric spheres of radii p and p + a. 

Note that: 
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i) in E n f(p, a) = 


-, so for fixed a > 0 and n, f(p, a) —> 1 as p 


ii) in E n , for fixed a > 0 and p. f(p, a) —> 0 as n —> oo; 

iii) in HP, for fixed a > 0 and n, f(p , a) —> e _ca as p —> cx), for some constant c > 0. 

To see why these phenomena interfere with a generalization to hyperbolic space 
of the method used earlier for Euclidean packings, consider the packings 7-p of the 
hyperbolic plane, by disks of fixed radius II. defined for each p >> 0 as follows. For 
each sufficiently large radius p » R, place disks of radius R on the circumference 
of a circle C p of radius p, so that: they cover all but perhaps one arc of the circum¬ 
ference; there are as many disks as possible without overlap; disks intersect only at 
points of the circumference. We now show that by taking R (and therefore p) large 
enough we can ensure that the fraction of the area of C p covered by the disks is as 
close to 1 as desired. 


The fraction of the area of C p which is in the annulus between C p and the 
concentric C p r for p > p' is of the order 1 — e p ~ p for large p. p r , and by taking 
0 << p — p' « R « p' « p we can ensure that most of this area is inside the 
disks of radius R - all except those regions outside pairs of touching disks of radius 
R and outside the circle C p >. plus the region near any uncovered arc of C p . But 
using the convexity of circles, the former regions are each contained in triangles of 
the form TUV (see Figure 7), so have negligible area, and another simple triangle 
argument applies to the region near any uncovered arc of C p . 

So by choosing R appropriately we could get almost all the area of C p to lie 
outside C'. 

Where in the Euclidean argument we used larger and larger cubes, in hyperbolic 
space we would use fundamental domains of cocompact subgroups of the isometry 
group Q of HP\ But we needed the fact, in Euclidean space, that the volume of 
the portion of a packing near the boundary of the fundamental domain would be 
negligible, while we see now that for large fundamental domains and large bodies, 
this is far from the case. In summary, where we used i) to show the existence of 
optimal packings in Euclidean space, in hyperbolic space we have instead iii), which 
for large spheres is approximately ii). This is the intuitive reason why there has 
been difficulty defining optimally dense packings in hyperbolic space for so long. 


II. Some responses to the problem 

We have summarized above the arguments that, for at least some packings 
in hyperbolic space, there seems to be no reasonable notion of density. We now 
consider how one might proceed with an analysis of optimal density. 

Two avenues of response that come to mind are: to replace the essentially global 
definition of density with something more local; or to find a way to define density 
for those packings where it is reasonable, together with a convincing argument for 
excluding the others. 

a. Completely Saturated Packings 

A packing is called completely saturated [FeK, FKK, Kup] if it is not possible 
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to replace a finite number of bodies of the packing with a greater total volume of 
bodies and still remain a packing. Intuitively, we think of a completely saturated 
packing as one that is locally densest. In [FKK], it was proven that any convex body 
of Euclidean space admits a completely saturated packing (and more generally any 
body with the strict nested similarity property) (see also [Kup]). In [Bow], it is 
proven that completely saturated packings exist for all bodies j3 in either Euclidean 
of hyperbolic space. The argument given there extends easily to finite collections 
B. 

Example 5 (a completely saturated packing with low density). 

As pointed out in [FKK], a completely saturated packing of Euclidean space 
is a densest packing. This is not true in hyperbolic space. In this example, we 
construct a pair of bodies /A, /A in El 2 such that two completely saturated periodic 
packings by {/A, fa} exist that have different densities. The reason, as we will see, is 
due to the fact that the length of the boundary of a region in the hyperbolic plane is 
comparable to its area. Let fa be the tile shown in Figure 8. It is a regular octagon 
with all interior angles equal to 27r/8. Let fai be the unique periodic tiling by fa. 
Let fa 2 be the tile shown in Figure 9. It is formed from fa by adding “protrusions” 
to some edges and “indentations” to others. We will assume that these protrusions 
and indentations are made so that they fit together but are narrow enough so that 
there is a region of finite area C\ in each indentation that cannot be occupied by 
a nonoverlapping copy of fa 2 unless it is occupied by a protrusion. Also we assume 
that each protrusion fits into a unique indentation. 

fa 2 admits a unique periodic tiling T 2 . Let fa be equal to fa with a small hole 
removed from its interior. Let V be the obvious periodic packing by fa (i.e. the 
one that comes from T 2 by removing a small hole from the interior of each tile). 
Since fa admits a periodic tiling, it is clear that the optimal density of {fa, fa} is 
one. Just as clear, is the fact that the density of V is area(fa)/area(fa) < 1. We 
will show that V is completely saturated (if the hole in fa is small enough). 

It is a standard fact of hyperbolic geometry that there exists a constant C 2 > 0 
depending only on the symmetry group of fa (and the fact that fa contains a 
fundamental domain for this group) such that for all finite subtilings T of fa, 
\dfa'\ > C 2 \T\ (by \dfa'\ we mean the number of edges contained in exactly one 
tile of fa' and by \fa'\ we mean the number of tiles in fa'). Since the hole in the 
interior of fa 2 can be made as small as we like, we may assume that area(fa) > 
area(fa) — C1C2/2. 

Suppose for a contradiction that V is not completely saturated. Then there 
exists a finite subpacking V C V and another finite packing V" such that (V — 
V) U V" is a packing and area( V") > are fa V). We may assume without loss of 
generality that V fl V" = 0. 

We claim that the number of edges of V that have protrusions on them coming 
from bodies of V is at least \dfa'\/2. So let e be any edge on the boundary of V. 
Let e = eo, ei,.., e n be the sequence of edges defined by for 1 < i < n, e^\ and 
are on a body of V and ej + i is the “opposite side” of e* in the sense that if has 



a protrusion on it (relative to the body containing both e* and e* + i) then is 
its corresponding indentation and vice versa. This sequence is uniquely defined and 
ends in an edge e n on the boundary of V. It is easy to see that if eo corresponds 
to an indentation of V' (i.e. eo has an indentation on it coining from a body of V') 
then e n corresponds to a protrusion and vice versa. Thus the claim is proven. 

Note that it is not possible for any body of V" to fill completely any indentation 
on the boundary of V — V (in fact a region of area at least C\ is always unfilled). 
Hence the total area of V" is at most 

area(P") <|P'|areo(/?i) - (Ci/2) \dV'\ 

<\V , \[area(f 3 1 ) - CiC 2 /2] 

<\V' \area(j3 2 ) 

—area( V'). 

This contradicts the choice of V". So V is completely saturated. The moral is 
that, in hyperbolic space, locally densest does not imply globally densest. 

b. Controlling pathological packings 

We now discuss an approach to density specifically aimed at controlling those 
packings, such as the above example of Boroczky, which pose difficulty in computing 
a reliable density. Even though the methods are also applicable to Euclidean space, 
the interests of this article make it natural to specialize the discussion from now on 
to § = HP. 

The key idea is to use a pointwise ergodic theorem of Nevo ([Nev, Thm. 1] for 
dimension n > 3; [NeS, Thm. 3] for n > 2), the conclusion of which is the existence 
of limits of the type (1) in the intuitive definition of density. The fact that such 
theorems only prove existence of the limit “almost everywhere” is not a defect, it 
is a feature, necessitated by examples such as that of Boroczky. 

We begin by reproducing some notation and results from [BoR]. Let d{-. •) be 
the usual metric on S, and let O be a distinguished origin. We suppose given a 
finite collection B of bodies f3j in S. Let Eg be the space of all “relatively-dense” 
packings of § by congruent copies of the /3j, that is, packings V with the property 
that any congruent copy of a body in B intersects a body of V. On Eg we put the 
following metric, corresponding to uniform convergence on compact subsets of S: 

deiVi, V 2 ) = supjh(B k nV!,B k nV 2 )i (17) 

k> i k 

where B k denotes the closed ball of radius k centered at the origin, and for compact 
sets A and C we use the Hausdorff metric 


h{A, C ) 


max{sup inf d(a,c), sup inf d(a. c)}. 
aeA c £C ceC aeA 


( 18 ) 


It is not hard to see [RaW] that Eg is compact in this metric topology, and that 
the natural action: (g , V) G Q x E# —> di'B) G Eg of the isometry group Q of S on 
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Eg is (jointly) continuous. Let M.(B) be the family of Borel probability measures 
on Eg. We call a measure // G A4(B) “invariant” if for any Borel subset E C Eg 
and any g G G, g(gE) — g{E). Let A4j(B) be the subset of invariant measures 
and A4f(B) the convex extreme (“ergodic”) points of A4/(B), all in their weak* 
topology, in which M.(B) and M.i(B) are compact. 

We will study these ergodic measures as a substitute for studying individual 
packings. As we will see, for any ergodic measure // G Aij(B) there is a set of 
packings Z of full //-measure such that for each V G Z, the orbit of V is dense in 
the support of //. So studying //. is a lot like studying a packing in Z. We will make 
this relationship more clear in what follows but first some examples. 

Suppose V is a “periodic” packing, i.e. the symmetry group Tp of V is cocom¬ 
pact in Q. We will construct a measure ftp G J\4f(B) whose support is contained in 
the orbit 0(V) = {gV \ g G Gj C Eg of V. 0{V) is naturally homeomorphic to the 
(metrizable) space Q/Tp of left cosets by the homeomorphism qp : 0 (V) —> Q/Tp 
with qp(gV) — gTp. There is a natural probability measure on Q/Tp induced by 
Haar measure on Q by the projection map Tip : Q —> Q/Tp. (Aside from an overall 
normalization the measure on Q/Tp can be defined on sufficiently small open balls 
B C g/Tp as the Haar measure of any of the components of TT-p 1 (L> j.) Hence qp 
induces a probability measure ftp on (){V). This measure can then be extended to 
all of Eg in the following way: gp{E) — fip[E fl 0( V)\ for any Borel set E C Eg. 
We will use the term “periodic measure” to denote any measure in Mj(B ) associ¬ 
ated in this way with the orbit of a periodic packing. It is not hard to prove from 
the uniqueness of Haar measure on Q that there is only one probability measure, 
with support in the orbit of a periodic packing, which is invariant under Q. 

Next, we define the density of an invariant measure. After the definition, we 
will show how the density of an invariant measure relates to the density of packings 
in its support. 

For any p G EI n we define the real valued function F p on Eg as the indicator 
function of the set of all packings V such that p is contained in a body of V. (The 
latter condition will sometimes be expressed as p G V.) 

Definition 1. For any invariant measure //. G A4r(B). the “average density” /!(//) 
is defined as J Eg F p (y ) clg(y). 

Note: the average density Diji) is independent of the choice of p. because of the 
invariance of the measure, so p is not needed in the notation. For convenience we 
sometimes use p — O. 

If V fl is a random packing with distribution // then the above definition states 
that the density of // is the probability that the origin is contained in a body of V t , ; . 

For periodic packings V there is an obvious notion of density using a funda¬ 
mental domain of Tp. The above definition of density coincides with this intuitive 
notion for such special V. 

Proposition 1 [BoR], If V is a periodic packing, D(y x ) is the relative volume of 
any fundamental domain for Y x taken up by the bodies of x. 
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We need the following notation. As usual we let Q denote the group of orien¬ 
tation preserving isometries of hyperbolic n-space HP (for some fixed n > 2). Let 
7 r : > HP be the projection map g —» gO where O is some distinguished point in 
EI n . Then we let B r denote the inverse image under 7 r of the closed ball of radius 
r centered at O. Finally let A g denote a Haar measure on Q. normalized so that 
A g{B r ) is the volume of the r’-ball in HP. 

We will use the following special case of Theorem 3 in [NeS] to relate the density 
of an ergodic measure to the density of (almost every) packing in its support. 

Theorem 1 [Nevo], If Q acts continuously on a compact metric space X such that 
there is a Borel probability measure //. on X that is invariant and ergodic under this 
action, then for every function / G L P (X, p) (1 < p < 00 ) there is a set Z of full // 
measure such that for every z G Z, 

[ f dp — lim \ [ f(gz)d\g(g). (19) 

Jx r ^°° A g {B r ) Jd t 

Actually we will use the following extension of this result. 

Theorem 2. Under the same hypotheses as the above theorem, the set Z may be 
taken to be invariant under Q. 

We will prove this result in the next section. Applying Theorem 2 to the 
function F p and using the proof of Prop. 2 of [BoR] we get 

Theorem 3. If p G Mj(J3) then there exists a set of packings Z . of full //-measure, 
such that for all p G HP and all V G Z 


.. vol[PnBJr) 
Inn --- 

r -00 vol[B p (r)\ 


D(C). 


( 20 ) 


Note that this implies that the (closure of the orbit of the) stripe model has 
measure zero with respect to every invariant measure. We will give another expla¬ 
nation for this fact in a later section. 

From example 4 we concluded that it is not possible, in general, to compute 
the density of a hyperbolic packing using a certain tiling associated to it. In spite of 
this we will show that it is possible to compute the density of an invariant measure 
using an associated space of tilings. 

Let E be a (compact, invariant) space of packings of HP. Let // be a Isom + (HP) 
invariant measure on E. Suppose that 0 is a space of tilings (of HP) and that there 
is an equivariant map <f> : E —> 0. For example, 0 may be the space of Voronoi 
tilings [FeK] corresponding to E. For V G E such that the origin is contained in 
a tile of 4>{V), let t p (V) denote the tile of <p(V) containing the point p G HP. We 
claim that for any p: 



voi[p n r p (v )] 
voI[t p (V)] 


dp(V). 


( 21 ) 
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Define a function / : H n x I" x S -> I by f(p , </, V) = 1/voI[t p (V)\ if p and q 
are both in the tile t p {P) and p is in a body of V (otherwise f(p, q, V) = 0). Define 
a measure v on HP x HP by 

u{E x F) — / / / f(p,q,V)dvol(p)dvol(q)dp,(V). (22) 

J E J E J F 

Since // is invariant, it is easy to check that for all g G Isom + ( HP), u(gExgF) — 
u(E x F). The mass-transport principle [BeS] implies that u(E x HP) = z/(HP x F) 
for any measureable E C HP. But it can easily be checked that u(E x HP) = 

vol(E)D(n ) and z/(HP x E) — vol(E) J vol[P r\T p (V)]/vol[T p (V)]dp,(V) (for any 

p). This proves the claim. For emphasis, we repeat that when // is an invariant 
measure we can compute its density with respect to local structures such as the 
Voronoi tilings. If fi-p is a periodic measure and 0 is the space of tilings by a 
fundamental domain of T-p then there is a natural equivariant map from the orbit 
of V to 0. The above result then yields Proposition 1. 

We now define optimality through measures. 

Definition 2. D(B) = sup^^^g) F>(//) will be called the “optimal density for 
£>”, and any ergodic measure ft G Aij(B) will be called “optimally dense (for £>)” 
if D{ji) — D(B). We define “optimally dense packings” a little differently than in 
[BoR], We say that a packing V is optimally dense if there is an optimally dense 
measure // such that the orbit of V is dense in the support of // and for every p G HP 
D(g) is equal to the limit of the relative fraction of volume in expanding spheres 
centered at p taken up by bodies of V. 

One of the main results of [BoR] asserts the existence of optimally dense mea¬ 
sures. Also, it was proven that for every // G Mf(B) there exists a set Z of full 
//-measure such that for every V G Z. the orbit of V is dense in the support of / /. 
Using Theorem 3 this implies 

Theorem 4. For any finite collection B of bodies there exists an optimally dense 
measure // on Eg, and a subset of the support of //, of full //-measure, of optimally 
dense packings. 

Note: There may be many optimally dense measures for a given B. 

In [Bow] it is proven that the set of completely saturated packings has full- 
measure with respect to any optimally dense measure //. In other words, if F /( 
is a random packing with optimally dense distribution // then F /( is completely 
saturated almost surely. Example 5 shows that the converse is false. 

A major result of [BoR] was that the set of all radii r such that there exists an 
optimally dense periodic measure for the sphere of radius r (in hyperbolic space) is 
at most countable. Thus most optimally dense sphere packings are complicated. 

c. What Invariant Measures Avoid 

Some packings, such as the Boroczky example, do not have a well-defined 
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density. We claim this is “due to” the fact that the closure of the orbit of such a 
packing has measure zero with respect to every invariant measure p . In this section 
we prove this statement and show other examples of packings that are not “seen” 
by invariant measures. 

Example 6 (Penrose’s binary tilings) 

Perhaps the most relevant to the discussion in section I is the B consisting of 
the body f3 shown in Figure 3. (This is a minor variation on the tile in [Pen], and 
a special case of tiles in [MaM].) We know copies of this body can tile HI 2 , and 
since limits in Eg of tilings will again be tilings, if there were any invariant measure 
/r G A4[(K) with support in the orbit closure of such a tiling it would clearly have 
density 1. However we can see there is no such measure as follows. First consider the 
slightly simpler, and better known, example of the natural action of the isometry 
group Q of HI 2 (namely Q = PSL 2 (M)) on the boundary A of HI 2 , instead of its 
action on the set of tilings. Assume there is a measure p on A invariant under 
Q. Any hyperbolic element gy, G Q has 2 fixed points in A, pi , p 2) and moves all 
other points towards one and away from the other. From its invariance under gy,. 
p{{pi,P 2 }) = 1- Then considering that any elliptic element g e G Q has no fixed 
points in A, and p must also be invariant under g e , we get a contradiction. So there 
are no probability measures on A invariant under Q. Going back to our space ©5 
of tilings by our body (3, consider the function / from 0g to A, which takes each 
tiling to the point “pointed to” by the protrusion on each body in the tiling. / 
is obviously continuous. If there were a probability measure p on the space ©5 of 
tilings, invariant under the action of Q, we could define a corresponding measure pg 
on A by Pf(E ) = p (/ _1 [E]). Since no such pf exists, this proves no such p exists. 

Now assume the optimal density for £?, D(B ), is 1, with an optimal measure p. 
For each R > 0 consider the function on Eg 


Ir(V) = [ Fo(gV)d\g(g), (23) 

vol(B R ) Jp R 

which gives the relative area of the ball B R covered by the disks of V. From the 
invariance of of p, 


f R (V) dp(V) = / Fo(P) dp{V) = D(p) = 1 , 


'£b 


'£ S 


(24) 


so /r(P) = 1 for //-almost every V G Eg. Letting R run through the positive 
integers, and intersecting the sets of full measure we get for each such R, we see 
there is a set of packings of full measure which are tilings. Since the closure of a set 
of tilings can only contain tilings, and the support of p must be invariant under Q. 
that support is contained in the set of all tilings of (3. But we saw above that there 
can be no such measure as p. and this proves that D(B) / 1. (Using modifications 
of this example it can be shown [Bow] that for every e > 0 there exists a body f3 
that admits a tiling of HP and D((3) < e.) 
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The formalism above leads one to assert that the densest packings of the body 
/3 have density bounded away from 1, even though one can tile H 2 with copies of f3. 
The “reason” for this is that there are no invariant measures which can “see” the 
tilings; they are a set of measure zero for every invariant measure on the space of 
all packings by f3. We explore the consequences of this using some of the examples 
we discussed earlier. 

Consider again the tile f3 shown in Figure 3. Congruent copies of /3 can only tile 
the plane (up to an overall rigid motion) as in Figure 4 (in which the little bumps 
on the tiles are not shown.) Construct the tile (3 of Figure 5 out of three abutting 
copies of (3. Now drill a hole in (3 1 producing the body j3o, as shown in Figure 6. Note 
that the packings of the plane by (3q obtained in the obvious way from the tilings by 
(3, are precisely the complements of the disk packings of Boroczky discussed above. 
The point is, although it might seem reasonable to assign a density of 1 to the 
tiling of Figure 4, that would seem to imply a well defined density to the packing of 
Figure 2, which we know is misleading. In other words, the meaningfulness of the 
density of the tiling of Figure 4 is unstable under arbitrarily small perturbations 
(drilling arbitrarily small holes). Notice that when we drill these small holes we 
turn the tiling into a mere packing, forcing us to give up the “simplicity” of the 
tiling, as a global object with seemingly obvious density, and leaving us to find some 
meaningful way to assign a density to the resulting packing. As we will see below, 
the difficulty in assigning a density to a packing, for instance congruent copies of a 
single body (3, can derive from the complexity of the set of rigid motions of (3 that 
define the packing. And in this sense a tiling is no simpler; treating it as a global 
object with an “obvious” density simply avoids coming to grips with the essential 
nature of the assignment of density for packings. 

In other words, the phenomenon whereby the “optimal” density can be less 
(even far less) than 1 for a body which can tile space, can be understood as related 
to the instability of the meaningfulness of the density of the tilings under removal 
of small holes in the tiles. This suggests that even for tilings one needs to keep 
track of the individuality of the tiles. In this example that amounts to noting the 
various sets of congruences used in producing the tilings; in some sense those sets 
of congruences are too complicated to be analyzed through our density formalism. 

We have shown that if T is a tiling whose orbit closure in the space of packings 
factors onto the space at infinity of the hyperbolic plane, then there are no invariant 
measures on the orbit closure of T. All of our examples of strange behaviour in the 
hyperbolic plane have, so far, been constructed using this principle. Could this be 
the only way of constructing such examples? 

Example 7. 

The following example is a variant on a simple construction. First consider 
two congruent regular all-right-angles octagons in the plane. Label their edges in 
clockwise order ei, e 2 , es and e[, e' 2 , .., e' 8 . By identifying e*, with e' k for odd k (by 
orientation reversing homeomorphisms), we obtain a sphere X’ with four open disks 
removed. If we then identify pairs of boundary components of X\ the resulting 
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object is a genus two surface. The covering space S of the genus two surface 
corresponding to the commutator subgroup looks like the boundary of a regular 
neighborhood of the standard Cayley graph of the free group on two generators. In 
this way, we obtain a tiling of S by regular all-right octagons. 

For the variation, we wish to distinguish a boundary component of X'. We 
do this by modifying the edges of the all-right octagon so that e 2 (and e' 2 ) has a 
protrusion and e 4 ,e 6 ,e 8 (and e^, eg, eg) have indentations. We call this tile r (see 
figure 10). We identify eg with e' k for k odd as before to obtain X as in figure 11. 
We want to obtain a tiling of S by X. For this picture the standard Cayley graph 
of the free group on two generators. Draw arrows on each of the edges so that each 
vertex has exactly one outgoing arrow. For each vertex v. let X v be a copy of X. 
If there is an edge with arrow pointing from v to w, then we identify the boundary 
component of X v that has the protrusion with one of the boundary components of 
X w that has an indentation. In this way, we obtain a tiling of S by X (and thus by 
T>- 

Note that the free group F 2 on 2 generators {a, b} acts naturally and isomet- 
rically on S. This action lifts to an isometric action of H 2 via the covering map. 
Though not relevant to what follows, note that this lift is unique up to postcompo¬ 
sition by rigid motions of H 2 . 

If we start a walk in T$ from some initial tile and follow the protrusions we get 
“closer” to a point on the ideal boundary of S. It is not too hard to see that this 
point does not depend on the initial tile chosen but only on the tiling Ts . Therefore, 
there is a map from the space of tilings of S by X (defined similar to the same way 
Eg is defined) to the ideal boundary of S that commutes with the action of F 2 . 
Since the ideal boundary does not admit an invariant Borel probability measure 
(for practically the same reason that A does not admit an invariant measure), 
neither does the space of tilings of S by X. 

Suppose that there exists an invariant measure // whose support is contained 
in the orbit closure of T. Then this measure pushes forward via the covering map 
to a measure /ig on the space of tilings on S by r. This measure //g is invariant 
under the action of F 2 but this contradicts the previous paragraph. 

Now suppose that there is an equivariant map a> from the orbit closure (){T) 
of T in E r to A the boundary at infinity of the hyperbolic plane. Let p — 4>{T). 
Since <t> is equivariant, the stabilizer of T must be contained in the stabilizer of p. 
However, the stabilizer of T is noncyclic (since it contains an isomorphic copy of 
the fundamental group of S which is noncyclic). By the theory of fuchsian groups, 
the stabilizer of T does not fix any point at infinity. This contradiction shows that 
<f> cannot exist. 

III. Proof of Theorem 2 

We need the following well-known fact (see 7.1.1.2 in [AVS]): 

Lemma 2. There exist positive real constants Ci and c 2 (depending only on the 


15 



dimension n ) such that 


lim vol[B R ]e ClR — c 2 . 


R—>00 


Corollary 1. For r > 0, 


V0l[B R - r 
R —>00 vol[B R ] 


lim 


- e~ Cir . 


if h-.g 


is Borel, let A-(h), A + (h) : Q 
A_(h)(g) — lim inf 


A + (h)(g) — lim sup 


R >00 vol{B R ) 

1 


U {± 00 } be defined by 
^ h(g'g) dXg(g') 


rRoo" vol{B R ) 


'B F 


'Bf 


h(g'g ) dXg(g') 


(25) 

(26) 

(27) 

(28) 


Lemma 3. If h is any nonnegative Borel function on Q then A+(h) and A-(h ) are 
continuous. 

Proof. 

Let < 71 , g 2 G G be such that the distance between g\0 and g20 is r. Then 


A+(h)(gi) = lim sup 1 

R — >oo vol{B R ) 

1 


'B, 


— lim sup , 

R^oo Vol{B R ) JBng- 1 

— lim sup — j I 

R^oo Vol{B R ) Jn- 1 [B R ( gi O)] 

> lim sup ————- / 

R —>00 vol(B R ) 


Kg'9i) dXg(g') 
h{g') dXg(g') 

Hg')dXg(g') 
h(g') dXg(g') 


(29) 


= lim sup 


vol(B R - r 


'tt 1 [B R — r (g 2 0)} 

) 1 r 


R^oo Vol(B R ) Vol{B R _ r ) J 7V -i[B R _ r .(g 2 0)] 


h(g') dXg(g') 


—e~ Cir A + (h)(g 2 ) 


Since g± and (72 are arbitrary, A + (h) is continuous. The proof for A-(h ) is similar. |fi 

Proof of Theorem 2. Let Go be a countable dense subset of G ■ Let Zq be as in 
Theorem 1 . At first we assume that / is nonnegative. Let Zf = Pliego 9~ 1 Zq- 
Since Zf is a countable intersection of sets of measure 1 , / / ( Zf ) = 1 . By definition, 
for all z G Zf and for all g G Go, gz £ Zq. For x G X, define h x : G —> IP by 
h x (g ) = f(gx). For 2 G Zf and g E Go, we have 

A+{h z )(g ) = lim sup 1 

vol{B R ) 


h z (g'g) dX G (g') 


= lim sup ———— 

R —► oo VOl(B R ) 

= [ fdg. 

J x 


'Br 


’B r 


f(g'gz) dX G {g') 


(30) 
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The last equation holds since gz G Z 0 . By the previous lemma, A + (h z ) is 
continuous. So the above equations hold for all g G Q. Similarly, A_(h z )(g) — 
f x fdn for all g G Q and z G Zf. So 


lim 

R—>oo 


1 

vol(B R ) 



f(g'gz)d\ G (g') 



(31) 


for all g G Q and all z G Zf. So the set Z = U ge ggZf satisfies the conclusion of the 
theorem and completes the case when / is nonnegative. 

In general, we set / = /_(_ — /_ where /+ and /_ are nonnegative. By the above, 
there are invariant sets Z + and Z- of full // measure satisfying the conclusion of 
the theorem for f + and /_. The set fl Z_ is invariant, of full /j measure 

and for all z G Zf, we have 


/ f = U-f-dg 

>x ■! x 

= [ f+d/-i- [ 

J X J X 

= to -4^ f 

R^oo vol(B R ) Jb 


— lim 


1 


r^oo voI(Br) Jr 

i r 


f- dg 

/' f+{gz)\ G {g) 
f-(gz) Xo(g) 


= lim 


R— >oo vol(B R ) 


— lim 


1 


R^oo Vol(B R ) 


1 Bf 


IB, 


f+{gz) - f-(gz) Xc(g) 
f(gz) x G {g)- 


(32) 


This proves the theorem. 


IV. Summary 

In [BoR] the notion of “optimal density” was defined for packings in hyperbolic 
space EP, as above, through the use of probability measures, on a space of packings, 
invariant under the congruence group of HP. The notion of an optimally dense 
packing was also introduced in [BoR], but not very successfully. The justification 
for that term was not well connected to a limit (1); we could only show there that 
for a set of packings of full measure, the limit (1) existed relative to expanding 
spheres centered about any countable set of centers. One advance in this paper is 
an extension of Nevo’s ergodic theorem, allowing us to extend this proof of existence 
to all centers in HP, allowing a more natural notion of optimally dense packings. 

Perhaps more significantly, we have also contrasted our approach to density 
with earlier approaches, and compared some key examples, showing the significance 
of certain structural features of the space of packings to the existence of well defined 
densities. 


Acknowledgments. We are grateful to Russell Lyons for pointing out to us some 
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Figure 1. Boroczky’s packing of disks 
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Figure 2. Boroczky’s packing with two tiles in dark outline 
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